Ab Initio Molecular Dynamics Approach to Quantitative Description of Ion Pairing in Water

Ion pairing of lithium and fluoride in water is described quantitatively using ab initio molecular dynamics simulations. We design a reliable computational protocol for evaluating the ion-ion potential of mean force using density functional based simulation methods. By comparison to classical molecular dynamics with empirical force fields we establish the statistical error of the procedure. We also check the accuracy of the electronic structure description by comparison to experimental structural data and to higher level calculations for model systems. The present approach not only points to deficiencies in force field calculations of potentials of mean force for difficult cases of high charge density ions like the aqueous lithium fluoride pair, but also allows extracting electronic information, such as the amount of charge transfer to solvent and its dependence on the ion-ion distance

Peptide salt bridge stability: From gas phase via microhydration to bulk water simulations

The salt bridge formation and stability in the terminated lysine-glutamate dipep- tide is investigated in water clusters of increasing size up to the limit of bulk water. Proton transfer dynamics between the acidic and basic side chains is described by DFT-based Born-Oppenheimer molecular dynamics simulations. While the desol- vated peptide prefers to be in its neutral state, already the addition of a single water molecule can trigger proton transfer from the glutamate side chain to the lysine side chain, leading to a zwitterionic salt bridge state. Upon adding more water molecules we find that stabilization of the zwitterionic state critically depends on the number of hydrogen bonds between side chain termini, the water molecules, and the peptidic backbone. Employing classical molecular dynamics simulations for larger clusters, we observed that the salt bridge is weakened upon additional hydration. Consequently, long-lived solvent shared ion pairs are observed for about 30 water molecules while solvent separated ion pairs are found when at least 40 or more water molecules hy- drate the dipeptide. These results have implications for the formation and stability of salt bridges at partially dehydrated surfaces of aqueous proteins

Dynamics of Electron Localization in Warm vs. Cold Water Clusters

The process of electron localization on a cluster of 32 water molecules at 20, 50, and 300~K is unraveled using ab initio molecular dynamics simulations. In warm, liquid clusters, the excess electron relaxes from an initial diffuse and weakly bound structure to an equilibrated, strongly bound species within 1.5 ps. In contrast, on cold, glassy clusters the relaxation processes is not completed and the electron becomes trapped in a metastable surface state with an intermediate binding energy. These results question the validity of extrapolations of the properties of solvated electrons from cold clusters of increasing size to the liquid bulk

Developing machine-learned potentials to simultaneously capture the dynamics of excess protons and hydroxide ions in classical and path integral simulations

The transport of excess protons and hydroxide ions in water underlies numerous important chemical and biological processes. Accurately simulating the associated transport mechanisms ideally requires utilizing ab initio molecular dynamics simulations to model the bond breaking and formation involved in proton transfer and path-integral simulations to model the nuclear quantum effects relevant to light hydrogen atoms. These requirements result in a prohibitive computational cost, especially at the time and length scales needed to converge proton transport properties. Here, we present machine-learned potentials (MLPs) that can model both excess protons and hydroxide ions at the generalized gradient approximation and hybrid density functional theory levels of accuracy and use them to perform multiple nanoseconds of both classical and path-integral proton defect simulations at a fraction of the cost of the corresponding ab initio simulations. We show that the MLPs are able to reproduce ab initio trends and converge properties such as the diffusion coefficients of both excess protons and hydroxide ions. We use our multi-nanosecond simulations, which allow us to monitor large numbers of proton transfer events, to analyze the role of hypercoordination in the transport mechanism of the hydroxide ion and provide further evidence for the asymmetry in diffusion between excess protons and hydroxide ions

Machine learning potentials for complex aqueous made

Simulation techniques based on accurate and efficient representations of potential energy surfaces are urgently needed for the understanding of complex systems such as solid–liquid interfaces. Here we present a machine learning framework that enables the efficient development and validation of models for complex aqueous systems. Instead of trying to deliver a globally optimal machine learning potential, we propose to develop models applicable to specific thermodynamic state points in a simple and user-friendly process. After an initial ab initio simulation, a machine learning potential is constructed with minimum human effort through a data-driven active learning protocol. Such models can afterward be applied in exhaustive simulations to provide reliable answers for the scientific question at hand or to systematically explore the thermal performance of ab initio methods. We showcase this methodology on a diverse set of aqueous systems comprising bulk water with different ions in solution, water on a titanium dioxide surface, and water confined in nanotubes and between molybdenum disulfide sheets. Highlighting the accuracy of our approach with respect to the underlying ab initio reference, the resulting models are evaluated in detail with an automated validation protocol that includes structural and dynamical properties and the precision of the force prediction of the models. Finally, we demonstrate the capabilities of our approach for the description of water on the rutile titanium dioxide (110) surface to analyze the structure and mobility of water on this surface. Such machine learning models provide a straightforward and uncomplicated but accurate extension of simulation time and length scales for complex systems

Chasing charge localization and chemical reactivity following photoionization in liquid water

The ultrafast dynamics of the cationic hole formed in bulk liquid water following ionization is investigated by ab initio molecular dynamics simulations and an experimentally accessible signature is suggested that might be tracked by femtosecond pump-probe spectroscopy. This is one of the fastest fundamental processes occurring in radiation-induced chemistry in aqueous systems and biological tissue. However, unlike the excess electron formed in the same process, the nature and time evolution of the cationic hole has been hitherto little studied. Simulations show that an initially partially delocalized cationic hole localizes within similar to 30 fs after which proton transfer to a neighboring water molecule proceeds practically immediately, leading to the formation of the OH radical and the hydronium cation in a reaction which can be formally written as H(2)O(+) + H(2)O -> OH + H(3)O(+). The exact amount of initial spin delocalization is, however, somewhat method dependent, being realistically described by approximate density functional theory methods corrected for the self-interaction error. Localization, and then the evolving separation of spin and charge, changes the electronic structure of the radical center. This is manifested in the spectrum of electronic excitations which is calculated for the ensemble of ab initio molecular dynamics trajectories using a quantum mechanics/molecular mechanics (QM/MM) formalism applying the equation of motion coupled-clusters method to the radical core. A clear spectroscopic signature is predicted by the theoretical model: as the hole transforms into a hydroxyl radical, a transient electronic absorption in the visible shifts to the blue, growing toward the near ultraviolet. Experimental evidence for this primary radiation-induced process is sought using femtosecond photoionization of liquid water excited with two photons at 11 eV. Transient absorption measurements carried out with similar to 40 fs time resolution and broadband spectral probing across the near-UV and visible are presented and direct comparisons with the theoretical simulations are made. Within the sensitivity and time resolution of the current measurement, a matching spectral signature is not detected. This result is used to place an upper limit on the absorption strength and/or lifetime of the localized H(2)O((aq))(+) species. (C) 2011 American Institute of Physics. doi:10.1063/1.3664746

Antimycobacterial and Photosynthetic Electron Transport Inhibiting Activity of Ring-Substituted 4-Arylamino-7-Chloroquinolinium Chlorides

In this study, a series of twenty-five ring-substituted 4-arylamino-7-chloroquinolinium chlorides were prepared and characterized. The compounds were tested for their activity related to inhibition of photosynthetic electron transport (PET) in spinach (Spinacia oleracea L.) chloroplasts and also primary in vitro screening of the synthesized compounds was performed against mycobacterial species. 4-[(2-Bromophenyl)amino]-7-chloroquinolinium chloride showed high biological activity against M. marinum, M. kansasii, M. smegmatis and 7-chloro-4-[(2-methylphenyl)amino]quinolinium chloride demonstrated noteworthy biological activity against M. smegmatis and M. avium subsp. paratuberculosis. The most effective compounds demonstrated quite low toxicity (LD50 \u3e 20 μmol/L) against the human monocytic leukemia THP-1 cell line within preliminary in vitro cytotoxicity screening. The tested compounds were found to inhibit PET in photosystem II. The PET-inhibiting activity expressed by IC50 value of the most active compound 7-chloro-4-[(3-trifluoromethylphenyl)amino]quinolinium chloride was 27 μmol/L and PET-inhibiting activity of ortho-substituted compounds was significantly lower than this of meta- and para-substituted ones. The structure-activity relationships are discussed for all compound

Dynamic inundation simulation of storm water interaction between sewer system and overland flows

Copyright © 2002 Taylor & FrancisThis is the Author's Accepted Manuscript of an article published in the Journal of the Chinese Institute of Engineers (2002), available online at: http://www.tandfonline.com/10.1080/02533839.2002.9670691An improved urban inundation model, coupling a 2D non‐inertia overland flow model with a storm water management model, is adopted to simulate inundation in urban areas. The model computes, not only the overland runoff and the water overflow through manholes where surface runoff exceeds the capacity of storm sewers, but also the bidirectional flow interactions between sewers and overland runoff. The model was verified by a typhoon event in Nov. 2000, which resulted in serious inundation in the Mucha area of Taipei City. The result shows that the present model indeed improves simulation accuracy over the earlier model, and can be used to provide a more reliable flood mitigation design

Water Quality Monitoring: A comparative case study of municipal and Curtin Sarawak’s lake samples

In this study, particle size distribution and zeta potential of the suspended particles in municipal water and lake surface water of Curtin Sarawak's lake were compared and the samples were analysed using dynamic light scattering method. High concentration of suspended particles affects the water quality as well as suppresses the aquatic photosynthetic systems. A new approach has been carried out in the current work to determine the particle size distribution and zeta potential of the suspended particles present in the water samples. The results for the lake samples showed that the particle size ranges from 180nm to 1345nm and the zeta potential values ranges from -8.58 mV to -26.1 mV. High zeta potential value was observed in the surface water samples of Curtin Sarawak's lake compared to the municipal water. The zeta potential values represent that the suspended particles are stable and chances of agglomeration is lower in lake water samples. Moreover, the effects of physico-chemical parameters on zeta potential of the water samples were also discussed